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Abstract. 

We present a study of sensitivity of the R-parity violating SUSY searches with the upgraded D0 detector in 
Run II of the Fermilab Tevatron, within a SUGRA framework. We considered the lightest neutralino as an 
ON ■ LSP that decays into a lepton and slepton (R-parity violating decay), resulting in 2e+ > 4 jets or 2/x+ > 4 jets 

ON ■ final state. The analysis, based on scaling of the Run I results, shows that squarks and gluinos with masses up 

to about 0.6 TeV could be probed with 2 fb _1 of Run II data. This work has been done in the context of the 
BTMSSM Working Group of the Run II SUSY/Higgs Workshop at Fermilab. 
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I PHYSICS MOTIVATION 



Recent interest in R-parity violating (RPV) SUSY decay modes is motivated by the possible high-Q 2 event excess 
at HERA [1]. When interpretation of the excess through first-generation leptoquarks was excluded by the D0 [2] 
and CDF [3] experiments, it was suggested [4] that such an effect could be explained via the s-channel production 
of a charm or top squark decaying into the e + jet final state. Both the production and the decay vertices would 
thereby violate R-parity. Although more recent data has not confirmed the previous event excess, and despite the 
combined analysis showed that the anomalous events reported by the HI and ZEUS experiments were unlikely to 
originate from the production of a single s-channel narrow resonance [5], interest in RPV signatures has not abatted. 
The CDF and D0 Collaborations have recently performed searches for RPV SUSY [7,8], and have set new mass 
■ limits on the RPV SUSY particles. Both experiments focussed their searches on the A' couplings, as motivated by 
the high-Q 2 HERA event excess. The results of the D0 searches are extended to the Run 2 case and the expected 
• sensitivity to the RPV couplings is discussed. 

D 



II D0 SEARCH FOR RPV NEUTRALINO DECAYS 



The D0 search for RPV SUSY considered the case of neutralino LSP which decays into a lepton and two quarks due 
to a finite RPV A' coupling (see Fig. 1). Both the electron and muon decay channels were considered, corresponding 
to what commonly referred to as A^ and \' 2i j couplings, respectively. The corresponding final states contain either 
2e or 2/i and at least four accompanying jets. Unlike at HERA, this search is not sensitive to the value of the RPV 
coupling, as long as it is large enough so that the neutralino decays within the D0 detector. That corresponds to 
A' > 1CU 3 , which gives a lot of room, given current indirect constraints [6]. 




FIGURE 1. RPV decay of a neutralino LSP into a lepton and two quarks. 



We assume that the neutralino (LSP) pairs are produced in cascade decays of other supersymmetric particles and 
use all SUSY pair production mechanisms when generating signal events. 

Signal events were generated within the SUGRA framework with the following values of SUSY parameters: Aq = 0, 
/i < and tan (3 — 2 (the results are not sensitive to the value of Aq .) Center of mass energy of the colliding beams 
was taken to be 2 TeV. ISAJET [9] was used for event generation. The acceptance and resolution of the D0 detector 
were parametrized using the following resolutions: SE/E — 2% © 15%/VE [GeV] (electrons), S(l/p)/(l/p) = 
0.018 © 0.008(l/p) (muons), and SE/E = 3% © 80%/v^ [GeV] (jets) and found consistent with the full detector 
simulation based on GEANT [10]. 
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FIGURE 2. Points in the (mo,Jni/ 2 ) SUGRA parameter space used to generate RPV events in the ee + 4 jets channel. 

Figure 2 shows the points in the (mo,m 1 / 2 ) SUGRA parameter space where signal Monte Carlo events were 
generated for the electron channel. Similar points were studied for the muon-decay channel. 

Ill SELECTION CRITERIA FOR THE DIELECTRON CHANNEL 

A multijet trigger was used for the analysis of Run 1 data. It was found to be nearly 100% efficient for the typical 
RPV signal. Since Run 2 trigger list will include a similar trigger, we assume trigger efficiency of 100% and do not 
perform any trigger simulations for the Run 2. 

The following offline selections were used: 

• At least two good electrons, the leading one with E T {e) > 15 GeV and the other one with E T (e) > 10 GeV; 

• Rapidity range | rj |< 1.1 (central calorimeter), or 1.5 <| rj |< 2.5 (end calorimeters) for all the electrons; 

• Energy isolation for the electrons: the EM energy in the R=0.2 cone about the center of gravity of the EM 
cluster, subtracted from the total energy in R=0.4 cone, should not exceed 15% of the EM energy in the R-0.2 
cone. 

• At least four jets with Ex(j) > 15 GeV and | r\ |< 2.5; 

• The dielectron invariant mass (M ee ) should not be in the Z-mass interval, ie, | M ee — Mz |> 15 GeV/c 2 . 

In the present analysis we have dropped the requirement on Ht = J^Exie) + J2^tU) , but retained all other 
offline criteria that were used in the previous analysis of data from Run I [8] . 



IV SELECTION IN THE DIMUON CHANNEL 



The following event selection requirements were used for the muon decay channel: 

• Two muons, the leading one with px > 15 GeV, and the other one with px > 10 GeV. 

• Rapidity range \r]\ < 2.3 for both muons. 

• Energy isolation requirement for both muons, i.e. the calorimeter energy accompanying the muon in a (ry <j>) 
cone of 0.4 should be consistent with that from a minimum ionising particle. 

• At least four jets with E T (j) > 15 GeV and |?7| < 2.5; 

V SIGNAL EFFICIENCIES 

The number of signal events expected can be written as: (N) — C ■ a ■ e, where (N) is the expected number of 
events for luminosity £, a is the cross-section, and e is the overall efficiency. The efficiency e can be split into three 
terms: e = e t ri g ■ £kin ■ £id- Here e t ri g is the trigger efficiency for the events that pass the offline cuts ( assumed to be 
100%), ekin is the efficiency for offline criteria, which includes kinematic, fiducial and topological requirements, and 
£id is the electron/jet identification efficiency. 

The efficiency for identifying jets is very high (> 95%) and is expected to stay the same in Run 2. 

Electron identification efficiencies in Run 1 were 80 ± 7% in the central (|?7| < 1.1) and 71 ± 7% in the forward 
(1.5 < < 2.5) regions [8]. These efficiencies were calculated for electrons with Ex(e) > 25 GeV, It drops by about 
30% for electrons with E T (e) = 10 GeV. 

The muon identification efficiencies used in Run 1 were 62 ± 2% in the central (\r]\ < 1.0) and 24 ± 4% in the 
forward (1.0 < |r/| < 1.7) regions [11]. These were calculated for muons with px > 15 GeV. For muons with 
10 GeV < pt < 15 GeV the efficiencies were 80% smaller on average [12]. 

In the present analysis we have taken the overall particle identification efficiency to be 0.90 ± 0.09 in each channel, 
independent of lepton Ex, primarily due to the expectation of a better tracker and muon spectrometer for the 
upgraded D0 experiment. 

VI BACKGROUNDS 

The main backgrounds are expected to arise from Drell-Yan production in association with four or more jets, 
dilcpton top-quark events, and QCD multijet events. The latter is the dominant background for the electron channel 
(followed by the Drell-Yan background) . In the case of muons, the background is dominated by the Drell-Yan and 
top pair production. We used Monte Carlo to calculate background from the first two sources, and data to estimate 
background from QCD jets. 

Background for the Run 1 analysis was estimated to be 1.8 ± 0.2 ± 0.3 (with 1.27 ± 0.24 from QCD and 0.42 ± 
0.15 ±0.16 from the other processes) for ~ 100 pb^ 1 of data. To extrapolate this number to the data set from Run 2, 
we have simply multiplied it by the ratio of luminosities to obtain 36±4±6 events. However, it is expected that due 
to the central magnetic field in the upgraded D0 detector, the probability of jets to be misidentified as electrons will 
be reduced by a factor of ~ 2 in Run 2. We have therefore considered a second scenario with the smaller expected 
background of 15 ± 1.5 ± 1.5 events. 

For the muon channel, the expected background has been scaled directly from the Run 1 analysis. We expect 
10 ± 1 ± 1 background events in Run 2. 

VII RESULTS 

In order to obtain the sensitivity of Run 2 in to RPV decays, we calculated the efficiency for signal for all the 
mass points shown in Fig. 2. Typical efficiencies, the signal cross section in the ee + 4 jets channel, and the expected 
event yield in 2 of data, for several representative (mo, ^1/2) points, are given in Table VII. Similar numbers 
are obtained for the muon channel. 

We use these efficiencies to obtain exclusion limits in the (m , m 1 / 2 ) plane at 95% CL, assuming that no excess of 
events will be observed above the predicted background. The exclusion contours for the electron and muon channel 
are shown in Fig. 3 and 4, respectively. Numerical values of the limits are summarized in Table VII. 



TABLE 1. Efficiency x BR (%), signal cross section and the expected 
event yield in 2 fb^ 1 of data, at various (mo, m^) parameter space points. 



mo 


mi/2 


Efficiency x BR (%) 


Cross section 


(N) 


(GeV) 


(GeV) 




(Pb) 


(in 2 fb- 1 ) 


60 


235 


7.9 ± 1.1 


0.16 


25.2 ±3.4 


60 


245 


8.3 ± 1.1 


0.08 


12.8 ± 1.7 


60 


255 


8.3 ± 1.1 


0.06 


10.5 ± 1.4 


100 


220 


6.1 ±0.8 


0.10 


12.2 ± 1.7 


100 


230 


7.0 ± 1.0 


0.08 


11.3 ± 1.5 


180 


240 


7.0 ±0.9 


0.05 


7.1 ±1.0 


320 


240 


7.1 ±0.9 


0.05 


6.9 ± 1.0 



TABLE 2. Lower limits on the squark and gluino masses from Run 2. 







Lower limit on 


rriq 


Lower limit on m§ 


Limit when rriq — nig 






( For any m§ 


) 


( For any rriq ) 




Electrons 


Run 1 




252 GeV 




232 GeV 


283 GeV 


Run 2 (Scenario 


I) 


430 GeV 




490 GeV 


490 GeV 


Run 2 (Scenario 


II) 


520 GeV 




575 GeV 


585 GeV 


Muons 


Run 2 




560 GeV 




640 GeV 


665 GeV 



It's worth mentioning that our analysis provides a conservative estimate of the sensitivity achievable in Run 2, 
since no formal optimization of the signal vs. background has been performed. We expect that a formal optimization 
can improve the sensitivity in the mass reach by 15-20%. 
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FIGURE 3. Estimated exclusion contour for Run 2 in the (mo, Jni/ 2 ) plane for tanfi = 2, Aq = 0, n < 0, from the ee + 4 jets 
channel. Scenario I corresponds to a background of 36 ± 4 ± 6 events (direct scaling from Run 1); scenario II uses the 
background of 15 ± 1.5 ± 1.5 events (scaling, but with improvements in the detector taken into account). 
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FIGURE 4. Estimated exclusion contour for Run 2 in the (mo,m 1( / 2 ) plane for tan/3 = 2, Ao = 0, a < 0, from the aa + 4 jets 
channel for background of 10 ± 1.0 ± 1.0 (direct scaling from Run 1). 



